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The Boltzmann equation framework for inelastic Maxwell models is considered to determine the 
transport coefficients associated with the mass, momentum and heat fluxes of a granular binary 
mixture in spatially inhomogeneous states close to the simple shear flow. The Boltzmann equation 
is solved by means of a Chapman-Enskog-like expansion around the (local) shear flow distributions 
for each species that retain all the hydrodynamic orders in the shear rate. Due to the anisotropy 
induced by the shear flow, tensorial quantities are required to describe the transport processes instead 
of the conventional scalar coefficients. These tensors are given in terms of the solutions of a set of 
coupled equations, which can be analytically solved as functions of the shear rate a, the coefficients 
of restitution ars and the parameters of the mixture (masses, diameters and composition). Since 
the reference distribution functions apply for arbitrary values of the shear rate and are not 
restricted to weak dissipation, the corresponding generalized coefficients turn out to be nonlinear 
functions of both a and Ors- The dependence of the relevant elements of the three diffusion tensors 
on both the shear rate and dissipation is illustrated in the tracer limit case, the results showing that 
the deviation of the generalized transport coefficients from their forms for vanishing shear rates is in 
general significant. A comparison with the previous results obtained analytically for inelastic hard 
spheres by using Grad’s moment method is carried out showing a good agreement over a wide range 
of values for the coefficients of restitution. Finally, as an application of the theoretical expressions 
derived here for the transport coefficients, thermal diffusion segregation of an intruder immersed in 
a granular gas is also studied. 

PACS numbers: 05.20.Dd, 45.70.Mg, 51.10.+y 


I. INTRODUCTION 


Granular media under rapid flow conditions are 
amenable to a fruitful modelization through a gas of in¬ 
elastic hard spheres (IHS) [l|. In the simplest model, 
the grains are assumed to be smooth so that the inelas¬ 
ticity is characterized through a constant (positive) co¬ 
efficient of normal restitution a < 1 that only affects 
the translational degrees of freedom of the grains. The 
case a = 1 corresponds to elastic collisions. Due to the 
kinetic-energy dissipation in collisions, energy must be 
externally injected to the granular gas in order to main¬ 
tain it in rapid flow regime (fluid-like description). In 
some cases, the system is driven into the flow through a 
(linear) shear field (simple or uniform shear flow, USE) 
where a steady state is achieved when the energy dissi¬ 
pated by collisions is balanced by the energy supplied by 
shearing work. The study of the rheological properties 
in the steady USE has received consequential attention 
in the past years ii , especially in the case of monodis- 
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perse granular gases. 

The USE state is defined by a constant density n, a 
uniform granular temperature T, and a linear velocity 
prohle Ux = ay, where a is the constant shear rate. In 
the steady state, the system admits a non-Newtonian 
description 0,0 characterized by shear-rate dependent 
viscosity and normal stress differences. An interesting 
problem is the analysis of momentum and heat trans¬ 
port in spatially inhomogeneous states close to the USE. 
The physical situation is such that the granular gas is 
in a strongly sheared state that deviates from the USE 
conditions by small spatial gradients. The response of 
the system to these perturbations gives rise to additional 
contributions to the momentum and heat fluxes, which 
can be characterized by generalized shear-rate dependent 
transport coefficients. Due to the mathematical difficul¬ 
ties met in obtaining those coefficients from the Boltz¬ 
mann collision operator for IHS [I|, the inelastic version 
of the BGK model 0 was considered to determine the 
above generalized transport coefficients 0 , 0 . On the 
other hand, explicit expressions for these coefficients were 
derived by assuming particular perturbations where the 
steady state conditions of the USE apply 0 , 0 . This al¬ 
lowed us to perform a linear stability analysis of the hy¬ 
drodynamic equations with respect to the USE state 0 
to get the conditions for instability at long wavelengths. 
The results derived for IHS from the BGK model has 
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been then revisited by considering a mean field version 
of the hard sphere system where randomly chosen pairs 
of particles collide with a random impact direction. This 
assumption, which yields a Boltzmann collision operator 
with a collision rate independent of the relative veloc¬ 
ity of the two colliding particles, opens the possibility 
of obtaining exact results for granular gases in the con¬ 
text of the Boltzmann kinetic equation. The above in¬ 
teraction model is referred to as the inelastic Maxwell 
model (IMM) and it has been widely considered 

by physicists and mathematicians alike in the past few 
years to unveil in a clean way the role of dissipation in 
granular flows. In particular, the use of IMM allows us 
to analytically determine the set of generalized transport 
coefficients around the USF state for general unsteady 
conditions 0. 

All the above results refer to monocomponent granu¬ 
lar gases. However, a real granular system is generally 
characterized by some degree of polydispersity in den¬ 
sity and size (granular mixtures). Needless to say, the 
difficulties for obtaining explicit expressions of the trans¬ 
port coefficients increase considerably when one consid¬ 
ers multicomponent systems since not only the number 
of transport coefficients is larger than for a single gas but 
they are also functions of more parameters such as com¬ 
position, masses, sizes and different coefficients of resti¬ 
tution. In the case of states close to the homogeneous 
cooling state, explicit forms of the Navier-Stokes trans¬ 
port coefficients have been derived for IHS [l^ by con¬ 
sidering the so-called first Sonine approximation while 
exact expressions of these coefhcients have been also ob¬ 
tained for IMM [l^. In the case of far from equilibrium 
states, the results for granular mixtures are more scarce. 
In particular, the rheological properties (shear stress and 
normal stress differences) of inelastic Maxwell mixtures 
under USF has been explicitly determined in terms of the 
parameters of the mixture (concentration, masses, diam¬ 
eters and coefficients of restitution) [2^ . As in the 

case of monocomponent granular gases [T^, the use of 
the Boltzmann collision operator of IMM allows in prin¬ 
ciple to determine the transport properties in a strongly 
sheared granular mixture without introducing additional 
and sometimes uncontrolled approximations. In addi¬ 
tion, as has been mentioned in previous papers, the re¬ 
sults derived for inhomogeneous states from IMM com¬ 
pare well (especially in the case of low order moments) 
with those obtained from IHS [n m m, m, showing 
the reliability of IMM to assess the impact of collisional 
dissipation in granular flows. 

The goal of this paper is to study mass, momentum 
and heat transport in a strongly sheared binary mixture. 
In this case and taking the USF state as the reference 
one, the set of Boltzmann kinetic equations for the mix¬ 
ture is solved by means of a Chapman-Enskog-like ex¬ 
pansion 0 around the distributions of each species. 
Since the above distributions hold for arbitrary values of 
the shear rate [2lj, the different approximations in the 
Chapman-Enskog method retain all the hydrodynamic 


orders in a. Thus, the non-equilibrium problem ana¬ 
lyzed here accounts for two kinds of spatial gradients: 
small gradients due to the (slight) perturbations to the 
USF and arbitrarily large shear rates due to the reference 
shear flow state. In this paper, we will restrict our cal¬ 
culations to first order (Navier-Stokes-like hydrodynamic 
order) in the spatial gradients of concentration, temper¬ 
ature and flow velocity. It is important to remark that 
although the form of the zeroth-order distributions is 
not known, we only need their second- and fourth-degree 
velocity moments to evaluate transport around USF. The 
use of IMM instead of IHS allows us to exactly get these 
moments without the explicit knowledge of This 

is perhaps the main advantage of considering Maxwell 
models (both elastic and inelastic). 

In the first order of the expansion, the mass flux is 
characterized by the second-rank tensors Dij (diffusion 
tensor), Dp^ij (pressure diffusion tensor) and DT,ij (ther¬ 
mal diffusion tensor), the pressure tensor is defined in 
terms of the fourth-rank viscosity tensor rjijki while the 
heat flux is given in terms of the second-rank tensors 
D'lj (Dufour tensor), (pressure energy tensor) and A^ 
(thermal conductivity tensor). The set of the above gen¬ 
eralized transport coefficients are nonlinear functions of 
the shear rate, the concentration and the mechanical pa¬ 
rameters of the mixture (masses, sizes and coefficients of 
restitution). The determination of the equations defin¬ 
ing these transport coefficients is perhaps the main goal 
of the present contribution. 

As in previous papers pertaining to IMM [2l|, MM, 
the velocity moments of the Boltzmann collision opera¬ 
tor are given in terms of a collision frequency vq- This 
parameter can be seen as a free parameter of the model 
that can be chosen to optimize the agreement with the 
properties of interest of the original Boltzmann equation 
for IHS. Thus, in order to correctly describe the velocity 
dependence of the original IHS collision rate, one usually 
assumes that the IMM collision rate is proportional to 
with /3 = 1/2. Here, we take jd as a, generalized exponent 
so that different values of /3 can be used to mimic different 
interaction potentials. We assume that t'o oc nT^, with 
/3 > 0. In the case /3 = 0, vq is independent of tempera¬ 
ture (model A) while when jd ^ Q, vq \s a monotonically 
increasing function of temperature (model B). Model A 
is closer to the original model of Maxwell gases for elastic 
collisions [2^, [2^ while model B with fd = 1/2 is closer to 
IHS. The possibility of having a general temperature de¬ 
pendence of vq{T) for inelastic repulsive models has been 
also introduced in the granular literature [l^, [ 13 , [13 ■ 
of the main features of model A is that the reduced shear 
rate a* = ajvQ (which is the relevant parameter measur¬ 
ing the departure from the homogeneous cooling state) 
does not change in time and so, a non-Newtonian hydro- 
dynamic regime (where a* and the coefficients of restitu¬ 
tion ars are independent parameters) is achieved for long 
times. In this regime, the combined effect of both control 
parameters on the (scaled) transport coefficients can be 
studied analytically for model A. This is a bonus feature 
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of this model that contrasts with the results derived for 
model B where only analytical results can be obtained in 
the steady state limit (namely, when viscous heating and 
energy lost by collisions cancel each other). 

Our results show that in general the generalized trans¬ 
port coefficients associated with the mass, momentum 
and heat fluxes are given in terms of the solutions of a 
set of coupled nonlinear differential equations. In the 
case of model B (/3 7 ^ 0), these equations must be numer¬ 
ically solved with the appropriate boundary conditions 
to obtain their hydrodynamic forms. On the other hand, 
for model A (/? = 0), they reduce to a set of coupled 
algebraic equations that can be analytically solved. This 
allows us to provide explicit expressions of these coeffi¬ 
cients in terms of the shear rate and the parameters of 
the mixture. This achievement is perhaps one of the most 
relevant results of the present paper since it extends to 
binary mixtures previous results obtained for monocom¬ 
ponent gases m Nevertheless, the above expressions 
still involve quite a tedious algebra due essentially to 
the intricate dependence of the velocity moments of the 
zeroth-order distributions fr^'^ on both the concentration 
and the (reduced) shear rate. Thus, an exploration of 
the full parameter space, in principle straightforward, is 
beyond the scope of this presentation and a reduced prob¬ 
lem will be addressed. Indeed, the tracer limit (namely, 
a binary mixture where the concentration of one of the 
species is negligible) is specifically considered to analyze 
the behavior of the diffusion coefficients. 

Some results for diffusion around USF have been previ¬ 
ously reported by one of the authors of the present paper. 
Thus, in the tracer limit and for steady state conditions, 
explicit expressions of the diffusion tensors Dij , Dp^ij and 
DT,ij were derived for IHS [H, [ 2 ^ by using Grad’s mo¬ 
ment method [s^. In the case of IMM, the tracer diffu¬ 
sion tensor Dij has been also obtained [ 2 ^ in the steady 
state where the (reduced) shear rate a* is coupled to the 
coefhdents of restitution The results derived here 
for IMM extend to finite concentration the previous at¬ 
tempts made for the diffusion tensors in the tracer limit. 
In addition, given that all the previous works [13, [H, HI] 
have been restricted to the steady state, the results re¬ 
ported in this paper for model A (see subsection IV Al) 
offers the possibility of assessing independently the influ¬ 
ence of both a* and a^s on the diffusion of intruders in 
a sheared granular gas. 

Finally, as an interesting application of the general re¬ 
sults, a segregation criterion based on the thermal diffu¬ 
sion factor is derived in the tracer limit. This criterion 
shows the transition between two different regions (up¬ 
wards and downwards segregation) by varying the differ¬ 
ent parameters of the system. This study complements 
a previous analysis carried out for IHS [^. Our results 
show that the form of the phase diagrams of segregation 
is quite similar to those obtained before for IHS. 

The plan of the paper is as follows. In Sec. m 
the Boltzmann equation for inelastic Maxwell mixtures 
is introduced and the USF problem is defined. The 


Chapman-Enskog-like expansion around the USF state 
is described in Sec. IHII while Sec. |TV] deals with the eval¬ 
uation of the generalized transport coefficients associated 
with the mass, momentum and heat fluxes. The tracer 
limit is considered in Sec. El to illustrate the dependence 
of the (scaled) transport coefficients on the reduced shear 
rate and the coefficients of restitution. In this limiting 
case, the diffusion coefficients are the relevant transport 
coefficients of the mixture. The dependence of the coeffi¬ 
cients Dij , Dpjj and Dxjj on both a* and is analyzed 
for model A for general unsteady conditions while steady 
state conditions are assumed to get the form of Dij for 
model B (/3 7 ^ 0) in order to compare with previous re¬ 
sults derived for IHS [l^j. Comparison shows in general 
good agreement even for strong dissipation. Thermal dif¬ 
fusion segregation is studied in Sec. I VII while a brief dis¬ 
cussion of the results reported in this paper is provided 
in Sec. IVHI 


II. INELASTIC MAXWELL MIXTURES UNDER 
SHEAR FLOW 


A. Inelastic Maxwell mixtures 


We consider a granular binary mixture modeled as an 
IMM. In the absence of external forces, the set of nonlin¬ 
ear Boltzmann equations for the one-particle distribution 
function /r(r,v, t) of species r (r = 1 , 2 ) reads 

(^ + V • V) /^(r,V;t) = ^ JrsMfrit),fs{t)] , (1) 

'' <5=1 


where the Boltzmann collision operator J^s [vij/rj fs] for 
IMM describing the scattering of pairs of particles is 

Jrs[^l\fr,fs] = [ dv2 [ d? (v^ )/, (v^) 

nsild J J 

-/r(vi)/^(v2)] . (2) 


Here, 


rir 


dv/^(v) 


( 3 ) 


is the number density of species r, is an effec¬ 
tive collision frequency for collisions of type r-s, = 
27r'^/Vr(d/2) is the total solid angle in d dimensions, 
and ars < I refers to the constant coefficient of restitu¬ 
tion for collisions between particles of species r with s. In 
addition, the primes on the velocities denote the initial 
values {'v'i^'v' 2 ] that lead to {vi,V 2 } following a binary 
collision: 


v'l = vi - jisr (1 -b (o’ • §12)0-, ( 4 a) 

V 2 = V2 -b firs (I + (o’ • gl 2 )o’ , 


(4b) 
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where gi 2 = vi—V 2 is the relative velocity of the colliding 
pair, S’ is a unit vector directed along the centers of the 
two colliding spheres, and ^.rs = 'mr/{mr + rris) where 
rur is the mass of a particle of species r. 

Apart from rir, the relevant quantities in a binary mix¬ 
ture at a hydrodynamic level are the flow velocity u and 
the granular temperature T. They are defined, respec¬ 
tively, as 

1 2 2 

u = - ^ PsUs = X! / dvm^v/s(v), (5) 

P s=l s=l 

2 2 

nT = '^ UsTs = ^ / dv^yVs(v). (6) 

S—1 

In Eqs. m and (ED, Pr = mrUr is the mass density of 
species r, n = ni + n 2 is the total number density, p = 
Pi + P2 is the total mass density, and V = v — u is 
the peculiar velocity. Equations (0 and ED also define 
the (mean) flow velocity and the partial temperature 
Tr of species r. The partial temperature measures 
the mean kinetic ener gy of species r. As confirmed by 
computer simulatio ns | 32 ll . experiments and kinetic 
theory calculations the global granular temperature 

T is in general different from the partial temperatures T^ 
(non-equipartition of energy). 

Eurthermore, the mass flux for species r is defined as 

jr =mr J dvV fr{v), ( 7 ) 

the total pressure tensor is given by 

2 f 

/ dvm.VV/,(v), (8) 

s=l d 

and the total heat flux is 

2 p 1 

q = d^r-msV'^Vfs{^r). ( 9 ) 

s=l 

In addition, the rate of energy dissipated due to collisions 
among all the species defines the cooling rate C as 

C = f dv . ( 10 ) 

r,s 

At a kinetic level, it is also convenient to introduce the 
partial cooling rates Cr, measuring the rate of energy lost 
by species r. They are defined as 

Cr = 'y ^ Cl'S = ~ rp ^ ' / dv TUrV Jrs[/r, /s], 
s ^ ^ s 

( 11 ) 

where the second identity defines the quantities Crs- The 
total cooling rate C is given by 

2 

C = '^ XslsCs, ( 12 ) 

s=l 


where Xr = Ur/n is the concentration (or mole fraction) 
of species r and 7 ^. = JC/T. 

As said in the Introduction, one of the main advantages 
of considering IMM is that the moments of the Boltz¬ 
mann collision operator Jrs[fr, fs] defined by Eq. Q can 
be exactly evaluated in terms of the distributions fr and 
A without the explicit knowledge of both distributions 
[25l |. This property has been exploited to determine the 
second-, third- and fourth-degree collisional moments for 
a monodisperse granular gas [35l |. In the case of mixtures, 
only the first-, second-degree and third-degree collisional 
moments [T^ have been obtained. Their explicit forms 
can be found in the above papers. 

The results obtained before apply regardless the spe¬ 
cific form of the effective collision frequencies ujrs- These 
frequencies are independent of velocity but depend on 
space an time through its dependence on density and 
temperature. On physical grounds, oc rig. As in pre¬ 
vious works on IMM [l3, [2l|, [ 2 ^, we will assume that 
UJrs OC Us TP, with ^ > 0. The case /3 = 0 (a collision 
frequency independent of temperature) will be referred 
as model A while the case (3^0 will be called model B. 
The collision frequencies cors can be seen as free param¬ 
eters in the model to optimize the agreement with some 
property of interest of IHS. Here, ojrs is chosen to get the 
same partial cooling rate (rs as for IHS (evaluated by 
using a Gaussian distribution for fr). With this choice, 
(jjrs can be written as 0,0 

f Cr-s f Sr + 0s\ ^ A(Q\ rr0 

U}rs=Xs{—j - ) ^ 0 ’ 1^0 = A{fj)nT^^, 

(13) 

where the value of the quantity A will be defined later 
(see subsection IVBl) . In Eq. (IT^ . = {<Jr + as)/2, and 

S = 1 


B. Uniform shear flow 

We assume that the mixture is under USE. This state 
is macroscopically characterized by constant partial den¬ 
sities, a uniform temperature, and a linear velocity profile 

u(?/) = ui(y) = U 2 (y) = ayx, (15) 

where a is the constant shear rate. This linear veloc¬ 
ity profile, in computer simulations, can convenientl y b e 
generated by the Lees-Edwards boundary conditions [36l| , 
which are simply periodic boundary conditions in the lo¬ 
cal Lagrange frame moving with the flow velocity [s^ . 
Since Ur and T are here uniform, then the mass and heat 
fluxes vanish and the transport of momentum (measured 
by the pressure tensor) is the relevant phenomenon. At a 
microscopic level, the USE is characterized by a velocity 
distribution function that becomes uniform in the local 
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Lagrangian frame moving with the flow velocity u, i.e., 
/s(r,v,t) = fs(y,t). In that case, Eq. ([T|) becomes [2^ 

^~aVyW,= MhJi] + Mhj2] (16) 

and a similar relation for / 2 . The relevant balance equa¬ 
tion in the USE state is the balance equation for the 
temperature. It can be obtained from Eq. m and its 
counterpart for species 2 ; it is given by 


tifies the lack of equipartition of the kinetic energy. Obvi¬ 
ously, 7 = 1 for any value of the shear rate and/or the co¬ 
efficients of restitution in the case of mechanically equiva¬ 
lent particles (mi = m 2 , cti = cr 2 , and an = 022 = 0 : 12 ). 
Beyond this limiting case, the temperature ratio clearly 
differs from 1. For model A, 7 is determined from the 
condition 


X2Pi%. + {d-l)Plyy 

P2,xx + (d — 1)^*2* j,j, 


( 20 ) 


f) On* 

= (17) 

where C* = C/^'o, a* = a/^'o, Pxy = Pxy/p, p = nT 
being the hydrostatic pressure. Equation (fT^ shows that 
the temperature changes in time due to the competition 
of two opposite mechanisms: viscous heating (shearing 
work) and energy dissipation in collisions. It is apparent 
that, except for model A (/3 = 0), the collision frequency 
vo{T) oc Td is an increasing function of temperature, and 
so a*{t) oc T[t)~d is a function of time. Consequently, 
the (reduced) pressure tensor P*y depends on time in 
the hydrodynamic regime only through its dependence on 
a*{t) [ 2 ^. Therefore, for /3 ^ 0, after a transient regime a 
steady state is achieved in the long time limit when both 
viscous heating and collisional cooling cancel each other 
and the mixture autonomously seeks the temperature at 
which the above balance occurs. In this steady state, 
the reduced shear rate and the coefficients of restitution 
are not independent parameters since they are related 
through the steady state condition 

a*P:y = -^C. (18) 

On the other hand, when /3 = 0, dta* = 0 so that the 
reduced shear rate remains in its initial value regardless 
of the values of the coefficients of restitution ars- As 
a consequence, there is no steady state (unless a* takes 
the specific value given by the condition (ITSl) ') and a* 
and ars are independent parameters in the USE problem. 
Moreover, it must be also noted that the results obtained 
in the steady simple shear flow state are universal in the 
sense that they apply both for model A and model B, 
regardless of the specific dependence of z/q on T. The 
rheological properties for a granular binary mixture of 
IMM in the steady state were obtained in Ref. [ 2 ^ , while 
a more detailed study on the rheological properties has 
been carried out in Ref. [2l[. In particular, the shear 
stress P*y can be written as 

P4 = (19) 

where rf is the (scaled) nonlinear shear viscosity of the 
granular mixture. The dependence of 7 * on both a* and 
ars has been thoroughly analyzed in Ref. [2lj for different 
systems (see for instance. Figs. 3, 4, 5 and 7 of [H). 

Apart from the rheological properties, an interesting 
quantity is the temperature ratio 7 = T 1 /T 2 , which quan¬ 


where = Pr^ij /p and 

Pr,^y = J dVmrV^VyfriV). (21) 

The expressions of the partial pressure tensors have 
been obtained analytically for model A in Ref. [2lj. For 
model B, the forms of must be determined after solv¬ 
ing numerically a nonlinear set of differential equations. 
On the other hand, analytical results for P^ can be ob¬ 
tained in the case of model B in the steady state (where 
models A and B yield the same results). In this situation, 
the temperature ratio is obtained by solving the equation 

[m 


X2QPlxy 

^iClPlxy' 


( 22 ) 


where C = Cs and 

Psr /H , \ dr “t" ds 

1- —{l + ars}— - . 

(23) 

Here, u}*g = Ursl^'o- When the expressions of C* and P/ 
are substituted into Eq. (l2^ . one gets a closed nonlin¬ 
ear equation for 7 whose numerical solution provides the 
dependence of the temperature ratio on the parameters 
of the problem. As expected, the extent of equiparti¬ 
tion violation is greater when the mass disparity is large. 
Moreover, the predictions of IMM for 7 compare very well 
(see for instance. Figs. 2 and 3 of Ref. 0 ) with Monte 
Carlo simulations for IHS for conditions of practical 
interest. This excellent agreement shows again the reli¬ 
ability of IMM to capture the main trends observed in 
sheared granular flows. 


C = 

Srs 


2oj: 


Psr (1 T O^rs) 


III. CHAPMAN-ENSKOG-LIKE EXPANSION 
AROUND USE 

Let us now perturb the USE by small spatial gradients. 
The response of the system to those perturbations gives 
rise to contributions to the mass, momentum and heat 
fluxes that can be characterized by generalized transport 
coefficients. Our objective is to determine the shear- 
rate dependence of these coefficients for inelastic Maxwell 
mixtures. 

In order to analyze this problem we have to start from 
the set of Boltzmann equations © with a general time 
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and space dependence. Let Uq = a • r be the flow veloc¬ 
ity of the undisturbed USF state, where the elements of 
the tensor a are Uij = aSixSjy. As expected [13, 
the disturbed state the true velocity u is in general dif¬ 
ferent from Uq, and hence u = Uq -I- du, du being a small 
perturbation to uq. As a consequence, the true peculiar 
velocity is now c = v — u = V — i5u, where V = v — Uq. 
In the Lagrangian frame moving with Uq, the Boltzmann 
equations 0 can be written as 


the local version of the shear flow distribution (refer¬ 
ence state). This type of Chapman-Enskog-like expan¬ 
sion has been already considered in the case of mono¬ 
component granular gases to get the set of shear-rate de¬ 
pendent transport coefficients of IHS 0, [1] and IMM [l3| ■ 
More technical details on this method can be found in the 
above references. 

In the context of the Chapman-Enskog method [40 |. 
we look for a normal solution of the form 




(24a) 


- oEy -I-(V -f Uo)-V /2 = J22[/2,/2] + >/2l[/2,/l], 

(24b) 

where here the derivative V/i- is taken at constant V. 
The macroscopic balance equations for the densities of 
mass, momentum and energy associated with this dis¬ 
turbed USF state are obtained from Eqs. (I24al) and (I24bl) 
with the result 


dtUr -I- uo • Vur + V • (nr<5u) = — ^ , (25) 

nir 

dtSui -I- QijSuj -I- (uo -I- (5u) • VSui = (26) 

^ridtT ^n(uo -I- du) • VT = -aP^y 

s—1 ^ 

(27) 

where the mass flux j^, the pressure tensor P, the heat 
flux q, and the cooling rate are defined by Eqs. 0, m, 
and m, respectively, with the replacement V —>■ c. 

We assume that the deviations from the USF state are 
small. This means that the spatial gradients of the hy¬ 
drodynamic fields are small. For elastic gases, the specific 
set of gradients contributing to each flux is restricted by 
fluid symmetry, Onsager relations, and the form of en¬ 
tropy production [s^. However, for granular gases, only 
fluid symmetry applies and so there is more flexibility in 
the representation of the heat and mass fluxes since they 
can be defined in a variety of equivalent ways depend¬ 
ing on the choice of hydrodynamic gradients used. Some 
care is thus required in comparing transport coefficients 
in different representations using different independent 
gradients for the driving forces. Here, the concentration 
xi, the pressure p, the temperature T, and the local flow 
velocity du are chosen as hydrodynamic fields. 

Since the system is strongly sheared, a solution to the 
set of Boltzmann equations (I24al) and (I24bl) can be ob¬ 
tained by means of a generalization of the conventional 
Chapman-Enskog method (40j | in which the velocity dis¬ 
tribution function of each species is expanded around 


/,(r,V,t) = /,[A(r,f),V], (28) 

where 

A(r, t) = {a:i(r, t),p{r, f), T(r, f), du(r, t)}. (29) 

This special solution expresses the fact that the space 
dependence of the reference shear flow is completely ab¬ 
sorbed in the relative velocity V and all other space and 
time dependence occurs entirely through a functional de¬ 
pendence on the fields A(r,t). The functional depen¬ 
dence (1^ can be made local by an expansion of the 
distribution functions fa in powers of the hydrodynamic 
gradients: 

fa [A(r, t, V] = /i°) (V) + /(i) (V) + • • • , (30) 

where the reference zeroth-order distribution function 
corresponds to the USF distribution function but tak¬ 
ing into account the local dependence of the concentra¬ 
tion, pressure and temperature and the change V —>■ 
V — (5u(r, t) = c. The successive approximations 
are of order k in the gradients of xi, p, T, and du but 
retain all the orders in the shear rate a. Here, only the 
first-order approximation will be analyzed. 

When the expansion (l30l) is substituted into the defi¬ 
nitions ©-(nni), one gets the corresponding expansions 
for the fluxes and the cooling rate: 

+ , P = P(0) + P(1)+... , (31a) 


q = q(0) + q(i) + ... , ( = ((«) + + .. . . (31b) 

Finally, as in the usual Chapman-Enskog method, the 
time derivative is also expanded as 

dt = ^ + • • •, (32) 

where the action of each operator d[^'’ is obtained from 
the hydrodynamic equations dsni-iETi). These results 
provide the basis for generating the Chapman-Enskog so¬ 
lution to the Boltzmann equations (I24ap and (I24bl) . 


A. Zeroth-order approximation 


Substituting the expansions ((5(1 - (I5^ into Eq. (|24aL 
the kinetic equation for f[^'^ is given by 

f(0) 


^( 0 ) AO) _ y ^ 
Ot h gy 


_ T rf(o) f(o)] I T rf(o) f(o)i 
— 'llll/l i/l J + 'll2[/l 1/2 1- 


(33) 
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To lowest order in the expansion the conservation laws 
yield 

= 0 , 

^ (34) 

5ui + aijSuj = 0. (35) 


Since is a normal solution, the time derivative in Eq. 
(1551) can be represented more usefully as 



where in the last step we have taken into account that 
depends on (5u only through the peculiar velocity c. 
Substituting Eq. (IMll into Eq. (1^ yields the following 
kinetic equation for 




T 


on 


( 0 )' 


dT 


— acy 


df 


(0) 


dcx 

(37) 


A similar equation holds for Note that Eq. dSil) 

and its corresponding counterpart for apply for both 
models A and B. The zeroth-order solution leads to 
= 0. Therefore, the most relevant velocity 
moments of the distributions are the partial pressure 
tensors and (defined by Eq. (1^ by replacing 
fr —>■ fr)- They can be obtained from Eq. (1571) and 
its counterpart for when one multiplies both equa¬ 
tions by rrirCC and integrates over c. The set of coupled 
differential equations defining the above partial pressure 
tensors have been derived in Ref. [2lj (see Eqs. (26)-(30) 
in this article). These equations can be explicitly solved 
for general unsteady conditions in the case of model A 
(/3 = 0) where the tensors Pj°^ and P^^^ can be expressed 
in terms of the parameter space of the problem (shear 


rate, coefficients of restitution, masses, diameters and 
composition). In the case of model B (/3 ^ 0), analytic 
results can be only obtained in the steady state (where 
we recall that the expressions of the partial pressure ten¬ 
sors are the same for both models). Beyond the steady 
state conditions, in order to have Pj°^ and P^^^ for model 
B one has to solve numerically the set of differential equa¬ 
tions obeying these partial pressure tensors. This task is 
beyond the objective of the present paper since we are 
mainly interested here in providing analytic results. 

As mentioned in the Introduction, the solution to Eq. 
dSZl) has not been obtained so far (even for model A 
where the collision frequencies ujrs are independent of 
the granular temperature) and hence, the precise form 
of the zeroth-order distribution is not known. How¬ 
ever, an indirect information on the behavior of is 
given through its velocity moments. As a matter of fact, 
only the second- and fourth-degree velocity moments of 
are required to determine the generalized transport 
coefficients associated with the first-order solution fr^\ 
As we will show in Sec. IlYl while the partial pressure ten¬ 
sors Pj°^ and P 2 °^ are involved in the evaluation of the 
diffusion coefficients Dp^ij and DT,ij and the vis¬ 
cosity tensor rjijki, the fourth-degree velocity moments 

^r%ke of (defined by Eq. (IBSI) ) are needed to get 
the coefficients and Xij associated with the heat 

flux. Although the forms of Pr°^ (r = 1, 2 ) are explicitly 
known for granular binary mixtures [2I| . the moments 


^r%ke only known for the special case of monodis- 
perse granular gases 2^. Given the difficulty of obtain¬ 
ing those moments from the true Boltzmann collision op¬ 
erator, one could consider a BGK-like kinetic model for 
granular mixtures to evaluate them. 


B. First-order approximation 

The analysis to first order in the gradients is worked 
out in Appendix!^ Only the final results are presented 
here. The distribution function is of the form 

/(i) = Ai ■ Vcci + • Vp + Cl • VT -f T>i : V(5u, (38) 

where the vectors Ml,Si,Cl}, and the tensor Pi are 
functions of the true peculiar velocity c. They are the 
solutions of the following set of linear integral equations: 


dp 


) {Pdp + Tdr) Ai- acy^^ + £Mi + M1A2 = A 


dcx 
2 a i9Pa 


(0) 

xy 


d dxi 


dxi 


Si 


/ 2aT dP, 


(0) 


xy 


\ dp dxi 


T 


dxi 


(39) 



















{pdp + TdT)B^- 


'2adP^f , ^(0) , 


d dp 


-p- 


dp 


- 




1 — ac. 


= Bi - 


( 0 )' 


( p(o) _ 

dp‘^ \ ^ dp 


^dc, 

p dp 


' £,\Bi + AA 1 B 2 


Cl, 


( 40 ) 


(pdp + TdT)Ci- 


P( 0 ) , HO) , M: ^ 

dp ^ dp dT dT 


C 1 CLCy 


dCi 

dcx 


pCiCi J^iC>2 — Cl H 


2 a 9Pa 


(0) 


xy 


d dT 


d((0)' 

' dT 


(41) 


— (P^p + TSt) Tim ~ O'Cy ~ aSkyTiM + CiVim T M-iT^m = Dim- 


(42) 


Here, Ai, Bi, Ci and are defined by Eqs. 

(IA6I) . respectively. Moreover, Ci and A4i are the lin¬ 
earized Boltzmann collision operators around the refer¬ 
ence USF state: 

CiX = - (Jii[/r, A] + Jii[A,/f)] + Ji2[A,/f ]) , 

(43a) 


Dtm = - 


Tmi 

P 


J dc Ci Cl j (c), 


(49) 


MiX = -Ji2[A°\x]. 


(43b) 


Pijk£ — ^ — 'klls I dc CiCjlds k^ir), 

s=l 

(50) 

D'l, = E = -S / dc c 2 c.A.i(c), 

(51) 


Tij — ^ ( Tgij, Tg ij — 2 ^ ^ Ci^s,j(c), (52) 


S=1 

2 


A similar equation for applies by setting \ M 2. 

It is important to note that for P = ^, Eqs. O-dllD 
are expected to have the same structure as that of the 
Boltzmann equation for IHS, except for the explicit form 
of the operators Lg and Mg. 

Once the form of the distributions is known, the 
first-order corrections to the mass flux j[^}, the pressure 
tensor and the heat flux can be obtained. They Upon writing Eqs. (j44|)-([53]) use has been made of the 
are given by 

.(1) mim2n dxi p ^ dp p ^ dT 


^ij — ^ ^ ^s,ijf ^s,ij — 2 / ^ 

s=l •' 


symmetry properties of Ar, Br, Cr and VrM- In general, 
the set of generalized transport coefficients defined above 
are nonlinear functions of the shear rate, the coefficients 
of restitution and the parameters of the mixture (masses, 


( 1 ) _ dSui 
drk ’ 


(44) sizes and concentration). 

(45) 


(46) 


where 


Dij — 


nm2 


JdcciAij{c), (47) 


IV. GENERALIZED TRANSPORT 
COEFFICIENTS 

This Section is devoted to the evaluation of the gen¬ 
eralized transport coefficients associated with the mass, 
momentum and heat fluxes. We consider each flux sepa¬ 
rately. 


A. Mass flux 


Dp,ij — 


pm I 


y dec* -Bij(c), 


(48) 


The constitutive form for the mass flux to first order in 
spatial gradients is given by Eq. (IT41) . To illustrate with 
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some detail the evaluation of the transport coefficients 
of the mass flux, let us consider the diffusion coefficients 
Dij, defined by Eq. gH). These coefficients can be ob¬ 
tained by multiplying both sides of Eq. (15^ by miCj and 
integrating over c. After some algebra, one arrives at 

{pdp + Tdr) Aj - yoD,, - a^kDkj 

p dxi dxi y d dxi ^ dxi J 

^ ( 2 aTdP^°'^ , P"DT ,^3 


mim2np \ dp dxi dxi j mim2nT 


In Eq. (l54l) use has been made of the results [I 

/" 1 , ^ A A A A \ mim 2 n 

J dc mid [L-iAi^j + MiA2,j) = - 


J dc midAij = 


Pi 


dP 


(0) 


dP 


(0) 


l,ij 


dxi 


dxi 


where 


P ^12 /, I N 

= —] —M 2 i (1 + ai 2 )- 
dp2 


(54) 

1 

D d^ij: 

(55) 

(56) 

(57) 


In the hydrodynamic regime, the diffusion tensor can 
be written as Aj = AAj where Dq = {pT/mim 2 i 2 o) 
and D*a is a dimensionless function of the reduced shear 

•'J 

rate a*, the coefficients of restitution ars, the mass ra¬ 
tio p = mi I m 2 , the ratio of diameters (Ti/(T 2 and the 
mole fraction xi. The dependence of D*^ on the pres¬ 
sure and temperature is through the reduced shear rate 
a* oc T^~^/p. Thus, 


{pdp + TdT)D,, = {pdp + TdT)DoD*^ 

= (1-/3)Aj-/3AA' 


dDi 

da* 


(58) 


Consequently, in dimensionless form, Eq. (1541) yields 
'2a 


-P* -hC 

XV ' ^ 


(1 - /3) A, - /3a* 


-a* D* = 

p OXl 

2a* dP. 


X D* — 


xy 


d dxi 


dP* 

dxi 

dC 


dD*^ 

da* 

2a* dP* 


d dxi 


- 12^0*1 


dC 

dxi 


D. 


T,ij ■ 


(59) 


The equations defining the (scaled) tensors D*^^j and 
-J can be obtained by following similar steps as those 
made before for D*y After some algebra, the results are 


‘^p* .. + g* 

d M A ^ da* 


2a* 


Ay a 


dC 


D* 

P^3 


^ A - Pl^ 


-fa* 


Pi dP* dPi, 
p da* 


da* 


- I - C* -f a*^^ 1 D* 

’ d da* ^ ^ ft-* ' 


/3 + Cj 

dP* 

„ xy 


da* 

dD 


(60) 


T,ij 


2a* 

~d 


X P4 + (l-/3)a* J +/3C + (l-/3) 




X A.,, + + a*AT,fe, = -(1 - PK ( 


. {PI 9p:j 


p da* 


da* 


+ (1 “ /3) 3 a 


d da* 


-C+a*^ DP 


da* 


(61) 


Upon writing Eqs. dSO]) and dSB) use has been made of 
the identities 


pdpP}°^ =p{p*-a*da^P*), 


TAAf =p(l-/3)a*A.A:-, 


pdpC^°^ -iyoa*da^C. 


= (/3 - 1)C(°) + (1 - P)i'oa*da^C- 


(62a) 

(62b) 

(63a) 

(63b) 


In the absence of shear field (a* = 0), P*j — 1 

PPj = Xrjr^ij, and Eqs. dSHl-dni]) have the solutions 
D*j = D*6ij, Dp^^j = DpSij and A.ij = D^^ij where 


D* = [12*^ - (1 - p)cr 


' d-fi dC (r,* , ri*\ 

(64) 


Here, (* = /vq, P*^j = P^fj/p, v*jj = vd/vq, a*j = 
iiijii 2 {), Dpj^j = Dp ij/DpQ, and Drp^^j = Dx^ij jDpQ where 
Dpo = ip/pi2o) and we recall that = aSixSjy. It must 
be noted that P*y and C* depend also on xi and a* 
through their dependence on the temperature ratio 7 = 
T 1 /T 2 . 


(65) 

D*T = {P-l)^D*p. (66) 

P'd 
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Equations (l64l) - (l66l) agree with the expressions derived 
for IMM (model B with /3 = i) in the Navier-Stokes 
hydrodynamic order M- Beyond the Navier-Stokes do¬ 
main (vanishing shear rates), in general Eqs. (I5^ - (ICT) 
are nonlinear differential equations that must be solved 
with the appropriate boundary conditions. However, in 
the case of model A (/3 = 0), Eqs. (|59ll - (|6T|) become a set 
of coupled algebraic equations that can be readily solved. 


B. Pressure tensor 


The pressure tensor is defined by Eq. (H5|) in terms of 
the coefficients Eq. (l50l) . These coefficients can 

be obtained from the integral equation (|42|) after multi¬ 
plying it by nirCiCj and integrating over velocity. The 
result is 


2 a 




P!^l^ + C^°n{pdp + TdT)m ,ijki 0>jpf]l,ipk£ 0'pk£]l,ijpi^ (j~llf]l,ijk£ '^12£]2,ijk£^ 




(0) 


dp 


2 

dp 


dp dT 


(67) 


The corresponding equation for P 2 ,ijki can be obtained The coefficients rjr^ijkt can be written as rjr^ijki = 
from Eq. (1571) by changing 10 2. Upon writing (157)) use {p/^o)''!* ijki- The dependence of on p and T is 

has been made of the result through a* so that 


dCTTLiCiCj iPyI — '^llPl^ijkH '^12‘k]2^ijki: 


where 


Til = 


Wii 


d{d -I- 2) 


(1 -|- aii){d -|- 1 — ail) 


( 68 ) (pdp + TdT)pr,ijki = (pdp + Tdr) —rj* ^ 

Vo ’ 

{1 o\ PP *^dr,ijki 

= (1 - P)T]r,ijke - a 


+ 2 —^// 2 i(l + ai 2 ) 
d 


^ M 2 i (1 + Q;i 2 ) 

dT 2 


, (69) 


Vo da* 

(71) 


a;i 2 Pi 2 n , \2 

—^ 2 i(l + “ 12 ) ■ 
d(d + 2) p2 


Thus, in dimensionless form, Eq. dSil) can be finally writ- 
(70) ten as 




(977* 


/ * * I 5 ( 5 * ** \ /** I * * \ 

\^ipVl,jpk£ ' ^jpVl,ipk£ ^pkVl,ijp£) \P~llVl,ijk£ "r ^12^2,ijk£) 

dp* 2 / dP* 

= {Pki - a*dlykd [Pl^J - /3a* 


da* 


(72) 


where t*- = Tij/t'o- 


dijkt^ where verifies the differential equation 


In the case of mechanically equivalent particles, 
Pli]l^^ = Plij!^-^ = Pij’ dlijke/xi = vlijke/^2 = 


2 

* 71* I 

^a Pxy + C 




jki 


(1 -/3)??rjfef -/3a 

/ * * . *: * :(! * \ 7y*/v-i* 

\^ipVjpk£ ' ^jpVipk£ ^pkVijp£) ^pVijk£ 
dP* 

= -a*^ 4 ^ - {S^kP;e + SjkP*,) 

2 / dP*\ 

+^{p:i-a*v;yki) [p:,-pa*^y 


(73) 
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where 


* (l + a)(d+l-a) 

=- WdTT) -■ 

Equation (1751) agrees with the results derived in Ref. 0 
for a sheared monocomponent granular gas of IMM. In 
the limit of vanishing shear rates (a* = 0 ), the solution 
to Eq. (1721) can be written as 


rihkt = 
where 


(75) 


Xi1i[t22 - (1 - /3)C*] - 2 : 272^12 
[tii - (1 - /3)C1[t22 - (1 - /3)C] - 'ri 2 T -21 ’ 


(76) 


It must be remarked that a non-equilibrium phase 
transition has been recently [i^ identified in the tracer 
limit for a granular binary mixture of IMM. This tran¬ 
sition refers to the existence of a region (coined as the 
ordered phase) where the contribution of tracer parti¬ 
cles to the total kinetic energy of the system is finite. 
However, the above (surprising) behavior has been only 
analytically found when the collision frequency ujrs is as¬ 
sumed to be independent of the temperature ratio (“plain 
vanilla Maxwell model”) and hence, it does not seem to 
exist for the more realistic version of the IMM considered 
here. The effects of the above transition on the Navier- 
Stokes transport coefficients of a granular binary mixture 
has been recently studied 


A. Model A 


* ^ 3^272^1*1 - (1 - /3)C1 - a:i7i7-2i , , 

^ [tii - (1 - I3)C][t22 - (1 - /3)C] - n 2 r 2 i' 

For model B with /3 = 1/2, Eqs. (TfSll - dTTl) are consis¬ 
tent with those previously obtained for the Navier-Stokes 
shear viscosity of an inelastic binary Maxwell mixture 
0 - On the other hand, except in the above two limit 
cases, Eq. (l72)l for 77* y and its counterpart for 772 ijkt. 
can be only solved analytically for model A (/3 = 0). 

The evaluation of the transport coefficients associated 
with the heat flux is more involved than the one carried 
out before for the mass and momentum fluxes. For the 
sake of brevity, only the final expressions of the differen¬ 
tial equations defining the coefficients D'fj, Lij and A^ 
are provided (see Appendix iBl) . 


V. TRACER LIMIT 


In the tracer limit, ~ and the relevant ele¬ 
ments of the partial pressure tensor P^*^^ admit simpli¬ 
fied forms (see Appendix IQ . In particular, in the tracer 

limit, 71 ~ 7 , = i9a;i7 = = da-C = 0 and 

= P^}j/xi. Taking into account these simplifi¬ 
cations, and for model A (/3 = 0), Eqs. (l60l) - (l6T]) become 


2 a* 


Pixy + C Df - = -xf^PP 


(78) 


2 a* 

~d' 


n* * ^P^xy 

^2,xy - a 


■2C* 


DP, - 

9Phj 


-<kD;,k, = ^wPl,^, - a* 

dPf,, (2a*^dPl,y 


H”Q- 


da* 


da* 


-C]D^ 


T,ij 1 


(79) 


The results obtained in the preceding Section apply for 
models A and B and give all the relevant information on 
the influence of shear flow on the mass, momentum and 
heat transport of a granular binary mixture. Accord¬ 
ing to these results, the set of generalized (dimension¬ 
less) transport coefficients {L>* , 11 * ^^, , 77 ^^,...} 

are nonlinear functions of the (reduced) shear rate, the 
concentration xi and the mechanical parameters of the 
mixture (mass and size ratios and coefficients of restitu¬ 
tion) without any restriction on their values. On the 
other hand, the evaluation of these coefficients (even in 
the case of model A where the results are analytic) is 
quite tedious due essentially to the complex dependence 
of the partial pressure tensors and the temperature 
ratio 7 on both the mole fraction xi and the (reduced) 
shear rate a*. Thus, for the sake of simplicity, we con¬ 
sider the tracer limit (a;i —>■ 0 ) where the mass flux is the 
relevant flux since the momentum and heat fluxes of the 
system (intruder plus gas particles) are the same as those 
previously obtained [l7j| for a monocomponent granular 
gas of IMM. 


2a* 

~d 


p* 

^2,xy 


dPjxy 

da* 


* 7~)* 

^T,ij ~ ^D^T,ij ~ ^ik-^T,kj 


= xifi a 


da* 


— a 


dPU, 

da* 


2a*2 dP 


2,xy 


da* 


-C D, 


(80) 


Here, C* = (1 — o\f)l2d and = (a;j[' 2 M 2 i(l + Q;i 2 ))/(i 
where ujp = <^ 12 !vq. Upon deriving Eq. (l78l) we have ne¬ 
glected the contributions coming from the tensors D* 
and ^ since both tensors are proportional to xi (and 
hence, they vanish in the tracer limit) while D*^ is in¬ 
dependent of Xi- In addition, the derivatives da*Pfij 
and dg^Pf,, appearing in Eqs. (TfSll - dM)) are obtained in 
Appendix [Dj 

As in the case of IHS 0 , the coefficients Dij decouple 
from the other ones and hence, they can be obtained 
straightforwardly. Their expressions are 


D* = 


'D 


2a* r)* _ /-^t: 

d ^ 2 ,xy S 





n* P* 
^ik^l,kj 

2a* p* 
d ^ 2 ,xy 
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FIG. 1: (Color online) Panel (a): Shear-rate dependence 
of the (dimensionless) coefficients D%^/Dq and D^yjDQ for 
d = 3, cri/(T 2 = 1, mxjmi = 2 and three different values of 
the (common) coefficient of restitntion a\ a = \ (solid line), 
OL = 0.8 (dashed red line) and a. = 0.6 (dash-dotted blue line). 
Panel (b): Shear-rate dependence of the (dimensionless) coef¬ 
ficients Dyy/D q and D*^/Dq for the same parameter values. 
These results pertain to model A (/3 = 0). 


In the steady state {{2a*/d)P 2 ^^y -(- (^* = 0), Eq. (1511) 
is consistent with previous results derived for IMM un¬ 
der shear flow [2^ . The remaining coefficients ^ and 
are coupled and they obey the set of simple al¬ 
gebraic equations dZHl)-®. As alluded to above, they 
are proportional to the concentration xi and thus van¬ 
ish in the tracer limit. Yet, it is of interest to nor¬ 
malize them by their vanishing shear rate counterparts 
(which are also proportional to xi), to study their depen¬ 
dence on parameters other than xi. In order to illustrate 
the shear-rate dependence of the set of transport coef¬ 
ficients Aij = {Hj* , H* jj , we consider a three- 

dimensional {d = 3) granular mixture. Also, to reduce 
the number of independent parameters, the simplest case 
of a common coefficient of restitution {a = 022 = 012 ) 
is studied. Thus, the parameter space is reduced to four 
quantities {(7i/a2,rni/m2,a,a*}. 

According to Eqs. (I78l) - (l80l) . we have that A^z = 
Azx = Ayz = Azy = 0 in agreement with the symme¬ 
try of the linear shear flow (fT5)l . Thus, there are five 
nonzero elements of the tensors A^: the three diago- 


FIG. 2: (Golor online) Same as Fig. [T] for D* yy/D* Q and 
Dp^zy/Dp^o (panel (a)), DT^yy/D^fi and DT^xy/Dxfi (panel 

(b)). 



FIG. 3: (Color online) Plot of the ratio Tj^yxylff as a function 
of the coefficient of restitution a in the steady state for a 
monodisperse granular gas (IMM with /3 = 0, model A). The 
solid line is the result for a three-dimensional system {d = 
3) while the dashed line corresponds to a two-dimensional 
system (d = 2). 


nal {Axx, Ayy, and Azz) and the two off-diagonal ele¬ 
ments {Axy and Ayx). The algebraic equations 
also show that the anisotropy induced by the shear flow 
yields the properties Axx ^ Ayy = Azz and Axy ^ Ayx- 
The equality P^^yy = Pi^zz implies Ayy = Azz- This is 
a consequence of the interaction model considered since 
Ayy ^ Azz for IHS 0. 
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The shear-rate dependence of the relevant elements of 
the diffusion tensors D *-, D* ■■ and -j has been plot¬ 
ted in Figs. [T] and [2] for d = 3, cri/(72 = 1, milm 2 = 2 
and three different values of the (common) coefficient of 
restitution. Here, the tensors have been reduced with re¬ 
spect to their values at zero shear rate, namely, D*j/Dq, 
Dp,ij/D*^o and g where 


D* — 

Uq — 


'D 


-c 


Dt, = (82) 


D;o = Xij (^1 - (zzB - 2C* + . (83) 

It can be seen that the influence of the shear flow on the 
diffusion coefficients is in general quite important. We 
also observe that the anisotropy of the system, as mea¬ 
sured by the difference D*^ — D*y grows with both the 
shear rate and collisional dissipation. As expected, the 
shear field induces cross effects in the diffusion of parti¬ 
cles. This is measured by the off-diagonal elements D*y 
Pyx): {D*y^) and D^^^y {D^^yJ. These coeffi¬ 

cients give the mass transport along the x (y) axis due 
to spatial gradients parallel to y (x) axis. All these co¬ 
efficients are negative in the region of parameter space 
explored. We see that, regardless of the value of a, 
the shapes of the off-diagonal elements are quite sim¬ 
ilar: there is a region of values of a* for which their 
magnitude increase with increasing shear rate, while the 
opposite happens for larger shear rates. With respect to 
the diagonal elements, they are monotonically decreasing 
functions of the shear rate (shear-thinning effect), except 
in the region of small shear rates. In addition, Figs. [T] 
and [H also show that, at a given value of a*, their values 
decrease with dissipation. 

It is also interesting to weigh the respective importance 
of the zeroth- and first-order contributions to the (non¬ 
linear) shear viscosity. The zeroth-order USF viscosity, 
ry*, is defined by Eq. (HI while its (dimensionless) first- 
order contribution is given by the coefficient rj^yxy 
the steady state and for mechanically equivalent parti¬ 
cles, the ratio V^yxy/v* can be obtained from Eq. (1751) 
for model A (/3 = 0): 


Vxyxy (1 T o) I 2As 

^ 2(d + 2) (z/;-C)(l+6A,)’ 


(84) 


where v* is defined by Eq. (1751) and A^ = ■ 

Figure 13] shows the ratio y^yxy/v* versus the coefficient 
of restitution a for spheres {d = 3) and disks {d = 2). 
For elastic collisions (a = I), a* = 0 in the steady state 
and so, rj^yxy = V*- Moreover, the zeroth-order solution 
rj* generically gives a significant contribution to the total 
non-Newtonian shear viscosity. Figure |3| also displays a 
pronounced shear thinning effect: when a decreases, the 
steady state gas departs more and more from equilibrium, 
and the resulting a* increases; this in turn leads to a 









FIG. 4: (Color online) Plot of the diagonal (dimensionless) co¬ 
efficients Dxx/Dq (panel (a)) and Dyy/Dg (panel (b)) as func¬ 
tions of the (common) coefficient of restitution a in the steady 
USF state for d = 3 in the cases axja^ = 1 and mxjm^ = 2 
(C) and = 2 and milm 2 = 4 (E). The solid fines cor¬ 

respond to the results derived here for IMM (models A and 
B) while the dashed lines are the results obtained for IHS 

[li,]!!. 


decrease of y^yxy/v*- Thus, the shear thinning effect is 
more marked for y*yxy than for y*. Indeed, Ref. M has 
shown that the USF viscosity itself, ry*, exhibits shear 
thinning. For further details dealing with the detailed 
behavior of y*, see Ref. [2l[. 


B. Model B: steady state conditions 

In model B the collision frequency vo{T) is an increas¬ 
ing function of temperature and hence, the (reduced) 
shear rate a*(T) = a/r'o(r) depends on time. Thus, in 
order to determine the diffusion coefficients one would 
have to solve numerically Eqs. (I59l)-(|6ll) in the tracer 
limit, discard the kinetic stageof the evolution and elimi¬ 
nate time in favor of a*{t) (j,l2^. An additional technical 
difficulty in the case of granular mixtures is that the dif¬ 
fusion coefficients depend also on the temperature ratio, 
that is itself time dependent through its dependence on 
a*{t). The integration of Eqs. (l59ll ~ (l6T]l is therefore a sig¬ 
nificantly more complex problem than for the monodis- 
perse system. On the other hand, given that the results 
derived for the rheological properties in a single granu¬ 
lar gas under USE indicate that the influence of the 
temperature dependence on vq on rheology is quite small. 
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FIG. 5: (Color online) Same as Fig. |4]for D%y/D q (panel (a)) 
and Dyj./DQ (panel (b)). 


one can consider the steady-state solution for model B, 
which is at any rate of interest in its own right. In this 
case, the condition (fTSl) applies and the solution to Eqs. 
(|5^ (IMI) can be obtained analytically in the tracer limit. 
Here, we focus our attention on the tracer diffusion ten¬ 
sor D*j whose expression in the steady state is universal 
since it applies for both models A and B, regardless the 
specihc dependence of vq on T. 

A previous comparison between IMM and IHS for this 
tensor was carried out in Ref. [^ . However, the (approx¬ 
imate) theoretical results for IHS considered in Fig. 7 of 
were obtained from a Grad’s solution [s^ where the 
weight distribution is Gaussian [2^ instead of the shear 
flow distribution (zeroth-order solution) [2^. In the com¬ 
parison ^rformed here, we will use the latter predictions 
of IHS which are expected to be more reliable than 
the other ones . In this sense, the present comparison 
complements to the one made before in Ref. [2^. 

As in previous works on IMM [l^, [ 2 O, and in order 
to compare the results between IMM and IHS, the pa¬ 
rameter A appearing in the definition of vq (see Eq. (TT^ ) 
is chosen as 


/ 2(mi -fTO2) 




mim2 


(85) 


With this choice, the partial cooling rates Cr (associ¬ 
ated with the partial temperatures T^) of IMM (with 
(3 = 1/2) are the same as those obtained for IHS (as 
evaluated in the Maxwellian approximation) (^ . The 
dependence of the set of tracer diffusion coefficients 


{D*y.,D*y,D*y,D*y.} on the (common) coefficient of 
restitution a is illustrated in Figs. 0] and [S] for two dif¬ 
ferent systems. We observe in general a good agree¬ 
ment between IMM and IHS, especially in the case of 
the coefficients D*^. and D*y. These coefficients measure 
mass transport in the flow direction (x axis). It must 
be pointed out that the discrepancies between both in¬ 
teraction models turn out to be more significant as the 
disparity of masses or sizes increases. 


VI. AN APPLICATION: SEGREGATION OF AN 
INTRUDER BY THERMAL DIFFUSION 

As an application of the previous results, this section 
is devoted to the study of thermal diffusion segregation 
of an intruder in a sheared granular dilute gas. Segrega¬ 
tion and mixing of dissimilar grains is one of the most 
interesting problems in granular mixtures, not only from 
a fundamental point of view but also from a more prac¬ 
tical perspective. This problem has spawned a number 
of important experimental, computational, and theoreti¬ 
cal works in the held of granular media, especially when 
the system is fluidized by vibrating walls [d^. In the 
case of sheared systems, some computational and exper¬ 
imental works in annular Couette cells [1^ have shown 
that granular materials segregate by particle size when 
subjected to shear. On the other hand, in spite of the 
relevance of the problem, much less is known on the theo¬ 
retical description of segregation in sheared granular sys¬ 
tems. Previous theoretical studies [d^ on the subject for 
dense systems have been based on a Ghapman-Enskog ex¬ 
pansion around Maxwellian distributions with the same 
temperature for each species. As mentioned before, the 
assumption of energy equipartition can be only justihed 
for nearly elastic gases which means small shear rates in 
the steady USE state. 

Thermal diffusion is caused by the relative motion of 
the components of a mixture due to the presence of a tem¬ 
perature gradient. As a result of this motion, a steady 
state is finally reached in which the separating effect aris¬ 
ing from thermal diffusion is balanced by the remixing 
effect of ordinary diffusion [d^. The new feature of our 
study is to assess the impact of shear flow on segrega¬ 
tion. Under these conditions, the so-called thermal dif¬ 
fusion factor A characterizes the amount of segregation 
parallel to the temperature gradient. However, due to 
the anisotropy induced by the shear field, a tensor A 
rather than a scalar A is needed to characterize segre¬ 
gation in the different directions. Here, for the sake of 
simplicity, we consider a situation where the tempera¬ 
ture gradient is orthogonal to the shear flow plane (i.e., 
dxT = dyT = 0, dzT ^ 0 and d^Uy = a = const.). In this 
case, the amount of segregation parallel to the thermal 
gradient is measured by the diffusion factor A 2 defined 
by the relation 


cHn T d In xi 
dz dz 


( 86 ) 
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reduces simply to dzPzz = 0. The pressure tensor has the 
form Pzz = pP*^{a*) and hence, the identity dzPzz = 0 
allows to express dzP in terms of dzT. The result is 

d\np {l-l3)a*ida,P:,)dlnT 

dz p;,-a*{d,.p*,) dz ■ ^ ’ 

Finally, the balance equation (1271) for the granular tem¬ 
perature yields 


(90) 


FIG. 6: (Color online) Phase diagram for segregation for a 
three-dimensional system (d = 3) and three different values 
of the (common) coefficient of restitution a = 0:22 = Q 12 : 
a = 0.9 (solid line), a — 0.8 (dashed line) and a = 0.7 (dotted 
line). 



FIG. 7: (Color online) Phase diagram for segregation for a 
three-dimensional system (d = 3) with 022 = 0.9 and ai 2 = 
0.7. The solid line corresponds to the theoretical prediction 
obtained from Eq. (1921) while th e sy mbols refer to computer 
simulations carried out in Ref. for IHS in the so-called 
LTu flow (Couette flow with uniform heat flux). 


If we assume that the bottom plate is hotter than the 
top plate {dzT < 0 ), then the intruder rises with respect 
to the gas particles if > 0 (i.e., dz Inxi > 0 ) while the 
intruder falls with respect to the gas particles if < 0 
(i.e., dz Inxi < 0 ). 

Our goal here is to determine A^ in a steady state 
with (5u = 0 and Si —>■ 0 (tracer limit) where the spatial 
gradients of T, p and Xi point in the z-direction. Under 
these conditions, the balance equation ( 1 ^ yields j^]. = 0 
where j^], is given by 


.( 1 ) „ dxi m 2 „ dp m2n2 ^ dT 

H,z- T P’^^dz T ^’^dz- 


(87) 


According to Eq. (1571) . the condition = 0 leads to 


dxi _ m 2 Dp^zz dp 7712712 Dt,zz dT 
dz miT Dzz dz miT Dzz dz 


( 88 ) 


In the steady state, the momentum balance equation (ES) 


Upon deriving (IM)) we have neglected the term dzQz since 
it is of second order in the gradients of xi, p and T. As 
said in section im Eq. (1901) establishes a relation between 
the (reduced) shear rate a* and the coefficient of restitu¬ 
tion 022 . 

Use of Eq. ((5^ into Eq. (1551) and substitution of Eq. 
(1551) into Eq. (1551) finally leads to 

^ Dt,zz - (1 - P)a* iP*z - a*ida.P*,))-^Dp,zz{da,P:,) 

' D*zz 

(91) 

where Dp^zz = x^ and Dt,zz = x-^ Equa¬ 

tion m provides the thermal diffusion factor A^ in terms 
of the diffusion coefficients D*.,^ and the (re¬ 

duced) pressure tensor P*^ and the derivative da* Pzz ■ To 
evaluate those quantities, we consider model A (/3 = 0) 
where P*^ and da*Pzz are given by Eqs. m and (inii, 
respectively. In addition, the explicit forms of the dif¬ 
fusion coefficients can be found by solving the set of al¬ 
gebraic equations (iTSI - dSOl) for 7 = j = z. The results 
clearly show that, while > 0 , the coefficients 
and do not have a definite sign. 

The condition A^ = 0 provides the segregation crite¬ 
rion for the upwards/downwards segregation transition. 
Thus, according to Eq. m and given that > 0 , the 
marginal segregation curve (A^ = 0 ) separating segre¬ 
gation towards the cold wall (A^ > 0 ) from segregation 
towards the hot wall (A^ < 0 ) is given by the condition 


{P:, - a*{da*P*zz)) D*t,zz = a*{da*P:z)D;,zz- (92) 

Although relation (IM)) holds for models A and B alike, 
the form of the phase diagrams for segregation (A^ = 0 ) 
depends on the interaction parameter /3, since the quan¬ 
tities da* Pzz j Pt zz Pp,zz differ in both models, even 

in the steady state. On the other hand, according to the 
previous results derived in the monodisperse case [ 13 , it 
is expected that the influence of P on segregation is very 
weak. 

Before analyzing the dependence of the parameter 
space on the form of the phase diagrams, it is instruc¬ 
tive to consider some limit situations. When the intruder 
and the particles of the gas are mechanically equivalent 
(toi = m 2 , CTi = 0-2 and 022 = 012 ), the two species do 
not segregate. This is consistent with Eq. dMl) since 
then = 0 so that = 0 for any value of 
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the (common) coefficient of restitution. Another inter¬ 
esting situation is the elastic limit (022 = 0^12 = 1; which 
implies a* = 0 in the steady state condition (IMll l. In 
this case, P*^ = 1 and = 0 so that, Eq. (1^ holds 

trivially for any value of the ratios mi/m 2 and o'i/(T 2 
(the intruder does not segregate). Beyond the above two 
limiting cases, the criterion (IM|) is rather complicated 
since it involves all the parameter space of the problem 
(mi/m2, cri/(T2, 022). 


Figure |6] shows the phase diagram in the 
{mi/m 2 ,cri/ 172 } plane for d = 3 and three differ¬ 
ent values of the (common) coefficient of restitution 
a 22 = Q!i 2 . For the sake of simplicity, we consider 
model A (/3 = 0). All zero contours of A^ pass through 
the point (1,1) since when mi = m 2 and cti = (72 all the 
species are indistinguishable for this system. We observe 
that when the intruder is smaller than the gas particles 
{aI < ( 72)1 the main effect of collisional dissipation 
(or equivalently the dimensionless shear rate a*) is to 
reduce the size of the down segregation region while the 
opposite happens when ui > <72- On the other hand, the 
impact of dissipation on the latter case is smaller than in 
the former case (when tri < 172 ) and the curves tend to 
collapse into a common one for sufficiently large values 
of the diameter ratio. It is also quite apparent that 
in general large intruders tend to move towards colder 
regions since the upwards segregation is dominant and 
occupies most of the parameter space. This conclusions 
contrasts with the results obtained for vibrated dense 
systems since intruders tend to move towards hotter 
regions as they get larger [i^. It is also important to 
remark that the conclusions drawn here for IMM agrees 
quite well with those obtained before for IBS (see Fig. 
5 of Ref. [HI), showing again the reliability of IMM to 
describe segregation in granular flows. 


As a complement of Fig. [51 Fig. [7] shows a phase dia¬ 
gram for a22 ^ Q;i2 (0122 = 0.9 and ai2 = 0 . 7 ) in the 
case /3 = 0. The theoretical results derived for IMM are 
compared here against recent computer simulations per¬ 
formed in Ref. in the so-called LTu state, namely, 
a steady state where the inelastic cooling is exactly bal¬ 
anced by viscous heating (as in the steady USF state) 
resulting in a uniform heat flux (sol . [Fij i . In the simula¬ 
tions, segregation is induced by a thermal gradient par¬ 
allel to the y-direction {dxT = d^T = 0 but dyT ^ 0) so 
that, the physical situation slightly differs from the one 
studied here theoretically. Nevertheless, when cti ~ 02 
the agreement with theory is good. More significant dis¬ 
crepancies appear when the intruder is larger than the 
gas particles since in this case the theory predicts that 
intruders only move towards hotter regions (upwards seg¬ 
regation) . This contrasts with simulation data since they 
still show a small region of downwards segregation. 


VII. CONCLUSIONS 

In conclusion, we have investigated the mass, momen¬ 
tum, and heat fluxes for a binary mixture of inelastic 
grains. The system is driven out of equilibrium by an 
imposed shear flow, which injects energy while dissipa¬ 
tive collisions between the grains act as an energy sink. 
A kinetic theory description was proposed, where the in¬ 
tractable Boltzmann equation is simplified in a Maxwell 
model fashion. Such models are in some cases sim¬ 
ple enough to be amenable to a full analytical solution, 
while remaining true to the key physical phenomena un¬ 
der scrutiny. In this respect, our model is not the sim¬ 
plest possible of the Maxwell family (the so-called “plain 
vanilla” approach), since the collision frequencies iOrs are 
taken to be the same as those found for IHS, see Eq. m- 
In this equation, a free parameter (5 is introduced. While 
/3 = 1/2 is the natural choice to reproduce inelastic hard 
sphere phenomenology, it also leads to a complex inter¬ 
play between shear and dissipation in the steady state. 
On the other hand, it is convenient to decouple these 
effects, which is possible when /3 = 0. We thus discrimi¬ 
nate two sub-models, referred to as model A and model 
B, having respectively j3 = 0 and ^ 0. Model A enjoys 
a larger parameter space than model B, which is at the 
root of the greater analytical tractability of the approach. 

Perturbing the USF, we analyzed the response of the 
fluid mixture, from which generalized transport coef¬ 
ficients can be identified. Due to the anisotropy in¬ 
duced by the shear, these quantities appear in tensorial, 
rather than scalar form. A Chapman-Enskog-like method 
around the shear flow distribution allows to derive the 
nonlinear differential equations obeying the set of gener¬ 
alized transport coefficients (see Sec. HVl) . Hopefully, in 
the case of model A (/3 = 0), the above equations become 
simple coupled algebraic equations whose solution unveil 
the dependence of transport coefficients on the key pa¬ 
rameters (shear, dissipation, concentration, size and mass 
ratio). 

To reduce the complexity of the problem and to illus¬ 
trate the impact of both shearing and collisional dissipa¬ 
tion on transport, we focussed on the limit where species 
1 has a much smaller concentration than species 2, the 
so-called tracer limit (xi 0). In doing so, mass trans¬ 
port becomes the relevant phenomenon to address, since 
momentum and heat fluxes coincide with their mono¬ 
component (inelastic) expressions. There are then in 
general 15 different diffusion transport coefficients, that 
couple the mass flux to the gradients of density, pres¬ 
sure and temperature. The simplified model worked out 
here reduces this number to 12 (two diagonal and two 
off-diagonal elements for each diffusion matrix). Our 
results hold for arbitrary values of the shear rate, and 
are not restricted to small dissipation. They show that 
shear driving notably affects mass transport. In addition, 
good agreement is reported between our Maxwell treat¬ 
ment and previously derived inelastic hard sphere results 
(here, the relevant view is that of model B, where in the 
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FIG. 8: (Color online) Plot of the scalar diffusion coefficient 
D*{a) (relative to its elastic value) and the zero shear-rate 
diffusion coefficient Dq (relative to its elastic value) as func¬ 
tions of the (common) coefficient of restitution ai 2 = 0^22 = ct 
in the steady state for d = 3 and two different systems: 
01102 = 2 and mi/m 2 = 2 (solid lines) and (ti/(J 2 = 1 and 
mi/m 2 = 3 (dashed lines). 


in agreement with the results obtained for IHS (see Fig. 
5 of Ref. [H). 

Finally, we wish to remark that, at the expense of a fur¬ 
ther simplification of the Maxwell model (addressing thus 
the aforementioned plain vanilla treatment [l^), it is of 
interest to study the impact on transport of a recently 
evidenced transition taking place in the intruder limit, 
where the minority species rather unexpectedly carries a 
finite fraction of the total system’s energy. Work along 
these lines is underway. 
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steady state, dissipation selects a unique reduced shear 
rate). Finally, we analyzed the segregation phenomenon 
of an intruder by thermal diffusion within our framework, 
which deciphers how shear impinges on the separating 
effect of a thermal gradient, opposed by the remixing ac¬ 
tion of diffusion. Our predictions are in fair agreement 
with inelastic hard sphere simulations for Couette flows 
sustaining a uniform heat flux. 

As pointed out in the Introduction, most of the works 
on granular mixtures dill have been derived by taking 
the so-called homogeneous cooling state as the reference 
state. In this case, the mass transport is characterized 
by the single scalar coefficients D, Dp and Dt (see Eqs. 
(IH^ (l66l) for IMM) instead of the tensorial quantities 
Dij^ Dp^ij and when the system is sheared. Al¬ 

though these scalar coefficients cannot be directly com¬ 
pared with the diffusion tensors obtained here, it would 
be interesting to gauge the effect of dissipation on dif¬ 
fusion in both situations (driven sheared case and freely 
cooling condition). In Fig. [8l we plot the scalar diffu¬ 
sion coefficient D* = \{Dl.^ -|- D*y 79*^) (which can 
be understood as a generalized diffusion coefficient in a 
sheared mixture) and the zero shear-rate diffusion coef¬ 
ficient Dq (defined in Eq. (1821) ) as functions of the (com¬ 
mon) coefficient of restitution in the steady state (where 
the results for these coefficients apply for models A and 
B) for xi —>■ 0 and d = 3. We have scaled both coef¬ 
ficients with respect to their elastic values. Given that 
the reference states in both descriptions (shear flow state 
against homogeneous cooling state) are quite different, 
there are significant quantitative differences between D* 
and Dq . On the other hand, the dependence of both co¬ 
efficients on dissipation is qualitatively similar since they 
increase as a decreases. This tendency is more impor¬ 
tant in the freely cooling case than in the sheared state, 


Appendix A: Chapman—Enskog-like expansion 


In this Appendix, some technical details on the deter¬ 
mination of the first-order approximation by means 
of the Ghapman-Enskog-like expansion are provided. In¬ 
serting the expansions (I5(I)) - (I5^ into Eq. (I24al) . one gets 
the kinetic equation for 




ai4 




- -f (V-t- uq) • vl/f \ 


(Al) 


The velocity dependence on the right-hand side of Eq. 
(EH) can be obtained from the macroscopic balance equa¬ 
tions (l^ - (l?7l) to first order in the gradients. Using these 
balance equations in Eq. (ini, one gets 


dfl 


( 1 ) 
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Di,ij{c) = p 
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dp 
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dT 
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dp 


(^Pij ^ O'Vxyij'j 


( 0 )' 


(A 6 ) 


Upon writing Eq. (IA 6 I) use has been made of the identity 
= 0 and the expression of the total pressure tensor 
Pj-p of the mixture 

‘■J 


(1) _ _ dSuk 

^ij - 'hjkt : 


(A7) 


where pijkt is the viscosity tensor. 

The solution to Eq. (IA2p has the form given by Eq. 
(I55)l . where the coefficients A-i, Bi, Ci, and Vi are func¬ 
tions of the peculiar velocity and the hydrodynamic fields 
xi, p, T, and (5u. The time derivative acting on these 
quantities can be evaluated with the replacement 


di°'>ViSuj = Vidi°^6uj = -ajk^Mk- (All) 

The corresponding integral equations O-dlll) can be 
obtained when one identifies coefficients of independent 
gradients in Eq. (IA2p and takes into account Eqs. (IA9I) - 
(lAllI) and the mathematical property 



where in the last step it has been taken into account that 
X depends on 5u through c = V — i5u. 


Moreover, there are contributions from acting on the 
pressure, temperature, and velocity gradients given by 


Appendix B: Heat flux transport coefficients 

The heat flux is defined by Eq. (H51) in terms of the co¬ 
efficients D'/j (Eq. d^), Lij (Eq. ([32])) and Ay (Eq. (1531)1. 
In order to determine them, we introduce the quantities 
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The generalized transport coefficients U"-, Ay and Ay 
are defined as 


2 

^ ^ ^s,kkij ■ 
s=l 

(B4) 

The differential equations verifying the (scaled) co- 
efficients + m- 2 )vQln, = 

Lr,zjki(jni+m 2 )vo/T and A*^^ = Xr,ijki{'rni+m 2 )vo/p 
can be obtained by following similar mathematical steps 
as those made for the other transport coefficients. The 
final results can be written as 
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In Eqs. (IB5ll - (|B7p . we have introduced the fourth-degree E — oj* ~ Q^ii) _ (^ _|_ 1 ) _ 

velocity moments of the zeroth-order distribution id{d + 2) pi 2 2d{d + 2) 


^l%ke = ^ y c,CjCkCifA{c), 
and the (dimensionless) quantities 
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tions for the coefficients ^jke A,ijki can be 

obtained from Eqs. (IB5l) - (jB7l) by changing 10 2. As in 
(BIO) the case of the previous transport coefficients, Eqs. (IB5I1 - 
(IB7I) become algebraic for model A (/3 = 0). Even for this 
model, the solution to the above equations requires the 
knowledge of the fourth degree moments whose 


8 d(d + 2){d + 4) d{d + 2)(d + 4) expressions are only known for a monodisperse granular 

x/r 2 i(l + 012 )^ [ 3 /r 2 i(l + 0112 ) - (d + 4)], (Bll) gas of IMM . 
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Appendix C: Rheological properties in the USF. 
Tracer limit 


The explicit forms of the (reduced) pressure tensors 
= ^ 2 ,y /?T- 2'72 and = Pi,ij/n 2 T 2 of the solvent 
(excess) and the solute (tracer) components, respectively, 
of a granular binary mixture (in the tracer limit xi —>• 0 ) 
of IMM under USF are provided in this Appendix. We 
consider here model A (/3 = 0) where the coefficients 
of restitution Urs and the (reduced) shear rate a* are 
decoupled. 

The non-zero elements of P^ij given by 


p* _ p* __ 

2,yy 2,zz ^ 2A(S) ’ 


1 + 2dA{a) 

1 -f 2A(S) ’ 

(Cl) 


where 


[1 + 2A(S)]" ’ 

2 (d + 2 ) a* 

(1 -I- 022)^ W22 ’ 


(C2) 


(C3) 


^22 = 



and A(a) is the real root of the cubic equation 


where 


F = 


^12 
d + 2 


d 21 (mi2 + M217) (1 + ai2f, 


(Cll) 


2llj* 

G = (1 + 0:12) [d -I- 2 - //2i(l + 0:12)] , (C12) 


2 a;* 

H = ^^^^ M21M12(1 + Ql 2 )^- (C13) 

Here, 

u;*2 = —— = a/mi 2 + M2i7i (C14) 


where 7 = T 1 /T 2 is the temperature ratio. The temper¬ 
ature ratio is determined from the constraint 


a;i7 


PI.. + {d- 1 )P*^^ 
d 


(C15) 


Since the collision frequency uj *2 is a nonlinear function of 
7, one then has to numerically solve Eq. (IC15|) to obtain 
the shear-rate dependence of the temperature ratio. 


Appendix D: Evaluation of the derivatives of the 
pressure tensors with respect to the shear rate. 
Tracer limit 


A (1 + 2 A )2 = ^, 
d 


namely 



'I , 

f 27.^0 M 

- cosh 


[6 

^ d )\ 


(C5) This Appendix addresses the evaluation of the deriva¬ 
tives da*P 2 ^ij and da*Pi^ij for model A (/3 = 0) needed 
to determine the tracer diffusion coefficients Py, Pp,y 
and Dt^i] in the tracer limit. In the case of the excess 

, . component, according to Eqs. (EH) and (IC2I) . one has 

[Gb) 


In addition, the long-time behavior of the granular tem¬ 
perature T{t) ~ T 2 {t) is T 2 (t) = T 2 ( 0 )e^“^^* where 

2 P 2 %a* C (1 + ^ 22 )^ ^^ l-ai 2 

d LO 22 ^22 d + 2 ‘2-duj22 

Upon obtaining the second identity in (IC7p use has been 
made of Eq. (IC2I) and the result C* = (1 — ct 22 )/‘^d. 

In the case of tracer particles, the relevant elements of 
P*^j can be written as [ 2 l| 


dP 


2,yy 


da* 


4A 

(I-f 2A)(l-f 6 A)’ 


(Dl) 


JPj.. 4(d-l)A 
“ da* (1-f 2A)(I-f 6 A)’ 


(D2) 


dPr 


2,xy 


da* 


1-2A 

(I + 2A)2(I + 6 A)“’ 


(D3) 




l,xy 


F + g^ 2 % - ^a*xT"Pi,.y 

^^22 + G 


Pi,.. = Xl 


(C 8 ) 


where use has been made of the identity 


a*^=2A 

da* 


1 + 2A 
1-f 6 A’ 


(D4) 


(C9) The calculations for the tracer particles are more intri¬ 
cate. Eirst, we derive both sides of Eq. (IU 8 I) with respect 
to a* to obtain the result 


(CIO) 


= a(o) -p 

da* yy yy da* ’ 


(D5) 
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where 

a(0) _ ^'2,yy) ~ Pl,yy^22i9a* 

yy ~ aw *2 + G 


(D6) 


The derivatives da-Pi ^v ^a^^i.xx can be also ob¬ 
tained from Eqs. (IC9|) and (ICIOI) . Their final forms can 
be written as 


= — 

yy Aw, 


Xi 


22 


G 


{F + HPlyy-Gx^^PlJ. (D7) 


^^ = AW+A«|1, (DIO) 

da* y y da* 


In Eqs. (ID6I) and (ID7E da*X = and we have 

introduced the quantities 


n = ——r + W]^ 2 ‘ 


2 wC| 


d + 2 


(D8) 


^.aS + a 2^, (Dll) 


_ /^21 ^ _ 1^21 jj 

Lt = ---TVCt, U = - - -TT-tl. 


2ujf ’ 


2u;f^ 


(D9) where 


A (0) ^iH{da*P2^xy) ^l,yy ^VV Pl,xy^22i^a- 

^y ~ Aw2*2 + G ’ 

^ x,H{da.PI,,) - 2Pl,y - 2a* A% - Pl,M 2 {da^ A) 

Aw|2 + G 

{f + HPlxx - 2 a*a:r'- Gx^^Plxx) ■ 


(D12) 


(D13) 


(DM) 


(D15) 


To close the problem, it still remains to get the quan¬ 
tity da*"f, which can be determined from the relation 
(IC15|) by taking the derivative with respect to a* in both 
sides of this identity. The result can be written as 


dj _ ^_i Ai°j + (d- 1)A^°^ 

9a* “ d- aL^ -{d- i)xr^ A^ ■ 


(D16) 
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